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neously broken time-reversal symmetry is studied in a model ceramic supercon- 
>- . ductor with anisotropic pairing symmetry. It is shown by Monte Carlo simulations 

that such a "chiral-glass" phase is truly stable even under the influence of screen- 



ing. Existence of frustration in zero external field, arising from the ci-wave pairing 
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Q^ \ symmetry of high-T c superconductors, is essential to realize this phase. Via a 



finite-size scaling analysis, critical exponents associated with the chiral-glass tran- 
sition are estimated to be v C g = 1-3 ± 0.2 and n C G = —0.2 ± 0.2. These values 
turn out to be close to those of the Ising spin glass. Phase diagram of the model 
is constructed and the implications to experiments are discussed. 
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1 Introduction 

Among a variety of macroscopic thermodynamic properties of superconductors, the 
type and the nature of possible thermodynamic phases is of central importance. 
For example, considerable attention has recently been paid to the possible phase 
of high-T c superconductors in applied magnetic fields. 1 For sufficiently random 
or dirty samples, the existence of a true thermodynamic phase with zero linear 
resistance was predicted (a vortex-glass phase). 2 Recent numerical works suggest, 
however, that the screening effects eventually destabilize the vortex-glass phase. 3 
In zero external field, by contrast, the only thermodynamic phase known to date 
either in clean or random superconductors is the standard Meissner phase. 

Meanwhile, recent experimental studies have revealed that cuprate high-T c su- 
perconductors have an anisotropic pairing symmetry, probably of the d x 2_ y 2-wave 
type. 4 Naturally, one may expect that such anisotropic nature of the supercon- 
ducting order parameter could give rise to novel thermodynamic properties not 
encountered in the conventional s-wave superconductors, possibly the appearance 
of a new thermodynamic phase. Unfortunately, this appears not to be the case 
in clean single crystals, although enhanced effects of thermal fluctuations give rise 
to a variety of interesting phenomena in high-T c superconductors. This is because 
the dj,2_ y 2-wave order parameter is characterized by a single phase variable of the 
condensate as in the conventional superconductors, and the resulting Ginzburg- 
Landau Hamiltonian, which is known to well describe various macroscopic prop- 
erties of superconductors, has essentially the same form as in the conventional su- 
perconductors. By contrast, in ceramic or granular high-T c samples, the situation 
may well differ because ceramic samples can be regarded as a random Josephson 
network and the anisotropic superconducting order parameter largely modifies the 
properties of the Josephson junction. 

One remarkable effect is the appearance of the V junction' characterized by 



the negative Josephson coupling, J < 0, across which the order parameter changes 
the phase by it. One important consequence of such ix junctions is the appear- 
ance of frustration even in zero external field. 5 Frustration is borne by the 'odd 
ring', a closed junction-loop having odd numbers of it junctions. On the basis of a 
single-loop model in which high-T c ceramics were regarded as an ensemble of nonin- 
teracting junction-loops, Sigrist and Rice 6 successfully explained the paramagnetic 
Meissner effect observed experimentally in certain high-T c ceramic samples. 7,8 Mag- 
netic moments spontaneously induced at the odd rings give rise to a paramagnetic 
response observed experimentally. 

It is also evident, however, that the paramagnetic Meissner effect itself is not 
directly related to a new thermodynamic phase, since an ensemble of noninteract- 
ing junction- loops is enough to cause the paramagnetic susceptibility, 6,9 just as an 
ensemble of free spins gives rise to the paramagnetic Curie law. In the present 
paper, we wish to address the question whether a new type of equilibrium phase 
closely related to the unconventional pairing symemtry of high-T c superconduc- 
tors is ever possible in certain ceramic samples. Indeed, one of the present authors 
(H.K.) recently proposed that such a novel thermodynamic phase might be realized 
in zero external field in certain ceramic high-T c superconductors. 10 This state is 
characterized by a spontaneously broken time-reversal symmetry with keeping the 
U(l) gauge symmetry, and is called a 'chiral-glass phase'. 11 The order parameter 
is then a 'chirality', which represents the direction of the local loop-supercurrent 
over grains. From a symmetry consideration, the nonlinear susceptibility was pre- 
dicted to diverge with a negative sign at the associated chiral-glass transition point. 
Frustration effect, which arises due to the random distribution of 7r junctions with 
the negative Josephson coupling, is essential to realize this phase. Note that in 
this chiral-glass state, unlike in the vortex-glass (gauge-glass) state under external 
fields, the phase of the condensate is not ordered, even randomly, on sufficient long 



length and time scales: The ordering occurs only in the loop-supercurrents, or in 
the chiralities. 

There are several experimental results which appear to corroborate the exis- 
tence of such a novel glassy zero-field phase in ceramic high-T c superconductors. 
Since the discovery of high-T c superconductors, it has been known that ceramic 
high-T c samples often exhibit a glassy behavior reminiscent of the spin glass. 12 ~ 14 
More recently, Leylekian, Ocio and Hammann observed via the noise and ac sus- 
ceptibility measurements that LSCO ceramic samples showed a glassy behavior 
reminiscent of the spin glass even in zero external field. 15 They also observed an 
intergranular cooperative phenomenon indicative of a glassy phase transition. It 
appears natural to interpret this cooperative phenomenon in terms of the pro- 
posed chiral-glass picture. More direct support of the chiral glass has recently 
been reported by the ac susceptibility measurements on YBCO ceramic samples. 
Thus, Matsuura et al observed a negatively divergent nonlinear susceptibility at an 
intergranular transition point, 16 consistent with the proposed chiral-glass picture. 

It should be emphasized here, however, that the theoretical analysis of Ref.10 
was based on an analogy to the XY spin glass, 17 and completely neglected the 
effects of screening (coupling of the condensate to fluctuating magnetic fields). 
Thus, the fate of the proposed chiral-glass phase in the presence of screening is 
not yet clear. It should be noted that the screening effect could be substantial in 
the intergranular ordering of ceramic high-T c materials, since the length unit to 
be compared with the penetration depth is the grain size (~ Ifim) rather than 
the short coherence length of the Cooper pair. As the screening effect makes 
the otherwise long-ranged interaction between the chirality (quenched-in half a 
vortex) short-ranged, one may wonder if it would eventually wash out a sharp 
phase transition and destabilize the chiral glass phase, just as it destabilizes the 
vortex-glass phase of type-II superconductors in a field. 



The purpose of the present paper is to investigate the stability of the hypo- 
thetical chiral-glass phase against this screening effect, and to determine whether 
the chiral-glass state really exists as a true thermodynamic phase or not. For 
this purpose, we perform extensive Monte Carlo simulations on a frustrated three- 
dimensional lattice model introduced by Dominguez et al, 18 in which the phase 
variables located at each site of the lattice are coupled to the fluctuating magnetic- 
field variables at each link. The model can also viewed as a spin-glass-type (i.e., 
random and frustrated) version of a U(l) lattice gauge theory. The same model was 
studied previously by Monte Carlo simulations with an interest in the behavior of 
linear 18 and nonlinear 19 susceptibilities. But, neither of these previous simulations 
was fully equilibrated, and thus, could not give any information whether a true 
equilibrium phase could exist or not. In the present paper, we perform equilibrium 
simulations on the same model based on an extended ensemble method recently 
proposed by Hukushima and Nemoto, 20 trying to determine the equilibrium prop- 
erties of the model. A part of the results was already reported in Ref.21. 

The remainder of this paper is organized as follows. In §2, we introduce our 
model and several physical quantities of interest. In §3, we give technical details of 
our Monte Carlo simulation. The results of Monte Carlo simulation are presented 
and analyzed in §4. By studying the Binder ratio associated with the chirality, we 
show that there indeed exists a stable chiral-glass phase with a spontaneously bro- 
ken time-reversal symmetry even in the presence of screening. Critical exponents 
characterizing the chiral-glass transition are determined with use of a finite-size 
scaling analysis. Phase diagram of the model is also constructed. Section 5 is de- 
voted to summary and discussion. We discuss in some detail the implications of the 
obtained results to the possible experimental detection of the chiral-glass state in 
ceramic high-T c superconductors. In particular, requirements for the appropriate 
samples as well as the method of detection are examined. 



2 Model 

We assume that weak links connecting the neighboring grains are distributed suffi- 
ciently dense, so that the system can be viewed as an infinite network of Josephson 
junctions which are not decomposed into finite clusters. Putting superconducting 
grains at the sites of a simple-cubic lattice, we model such ceramic superconductors 
by a three-dimensional lattice model of Josephson-junction array with finite self- 
inductance. Neglecting the charging effects of the grain, we consider the following 
zero-field Hamiltonian, 18 ' 19 

H = - £ J y 008(61, - 0j - A tJ ) + ^A 2 £(V x A)\ (2.1) 

<ij> ZJ ~ Z7r p 

where $i is the phase of the condensate of the grain at the i-th site of a simple 

cubic lattice, A being the fluctuating gauge potential at each link of the lattice, 

0o the flux quantum, Jy the Josephson coupling between the i-th and j'-th grains, 

and the lattice curl V x A is the directed sum of A^s round a plaquette. C is the 

self- inductance of a loop (an elementary plaquette), while the mutual inductance 

between different loops is neglected. The first sum is taken over all nearest-neighbor 

pairs, while the second sum is taken over all elementary plaquettes on the lattice. 

Fluctuating thermodynamic variables to be summed over are the phase variables, 

6i, at each site and the gauge variables, A t j, at each link. 

In this model, quenched randomness occurs only in the distribution of the 

Josephson coupling, Jy, which is assumed to be an independent random variable 

taking the values J or —J with equal probability (±J or binary distribution), 

each representing and 7r junctions. We also assume that Jy is independent of 

temperature and magnetic field. This assumption is more or less justified when 

the intergranular ordering occurs at a temperature considerably lower than the 



superconducting transition temperature of the grain. While our simulation is per- 
formed for this particular distribution of J^, one could expect from experience 
in spin-glass studies that the results would be rather insensitive to the details of 
the distribution, e.g., a slight asymmetry between ± J or the detailed form of the 
distribution. 

The bare Josephson penetration depth in units of lattice spacing is given by 



Ao = 1/VA (2.2) 

where C is the dimensionless inductance defined by 

I = (27T/(f)o) 2 JC. (2.3) 

Thus, larger inductance corresponds to smaller penetration depth with enhanced 
effects of screening. In the limit C — ► 0, or Ao — > oo, the screening effect becomes 
negligible and one recovers the XY spin-glass Hamiltonian of Ref.10. In the op- 
posite limit C — » oo, on the other hand, it can be shown that the model reduces 
to the noninteracting one. 19 Therefore, the system remains in the disordered state 
even at T = in this limit. 

The symmetry property of the Hamiltonian (2.1) was analyzed in detail in 
Ref.19. Contrary to the vortex-glass (gauge-glass) Hamiltonian, the Hamiltonian 
(2.1), defined in zero field, keeps the Z2 time-reversal symmetry in addition to 
the U(l) gauge symmetry. Frustration arises from the random distribution of 
both positive and negative Josephson couplings. This should be contrasted to 
the vortex-glass (gauge-glass) problem, where the associated Hamiltonian does 
not possess the time-reversal symmetry due to external magnetic fields, while the 
frustration arises from the magnetic field, not from the Jij. 

The Hamiltonian (2.1) can also be written in the dimensionless form as 

n = 7 = - E 4 oos(fc - 9, - A tj ) + -lE(Vx A) 2 , (2.4) 



where Jy = Jij/J is the dimensionless quenched random variable, taking values 
±1 with equal probability. 

The local chirality may be defined at each plaquette by the gauge-invariant 
quantity, 19 

k p = 2~ 3 / 2 ]T Jij sm(6i - 6j - Aij), (2.5) 

where the sum runs over a directed contour along the sides of the plaquette p. 
Physically, the chirality, k p , is a half (n) vortex, being proportional to the loop- 
supercurrent circulating round a plaquette p. If the plaquette p is frustrated, the 
local chirality k p tends to take a finite value, its sign representing either clockwise 
or counterclockwise orientation of circulating supercurrent. If, on the other hand, 
the plaquette is unfrustrated, it tends to take a value around zero. Note that 
the chirality is a pseudoscalar in the sense that it is invariant under global U(l) 
gauge transformation, 9i — > 9* + Ad, A^ — »■ A^, but changes its sign under global 
Z2 time-reversal transformation, 8i — >■ —$i, A^ — >■ —A^. Due to this symmetry 
property, chirality can be regarded as an order parameter of the chiral order. 

Induced local flux or magnetization threading a plaquette p is given in the 
dimensionless form by, 

f P = % % = V x A. (2.6) 

Flux is also a pseudoscalar like chirality, whose sign represents the direction of 
the induced magnetic moment threading the plaquette p. Total magnetization per 
plaquette along the z-axis is given by 

1 



m 



AnSN 



E %, (2-7) 



where S is the area of a plaquette and the sum is taken over all N p plaquettes on 
the < xy > plane of the lattice. The corresponding dimensionless quantity, rh, can 



be defined by 



™ - ~7T m = W S fp- ( 2i 

18,19 



The linear susceptibility, x> is given by 

dm n$N„ , 1 , 

dif £ L J 4vr v ; 

where if is an external magnetic field, /3 is the dimensionless inverse temperature 
defined by (3 = J/ksT, < ■ ■ ■ > represents a thermal average, and [■ • •] represents 
a configurational average over the bond distribution. The linear susceptibility is 
dimensionless in cgs units. The nonlinear susceptibility, xi, or its dimensionless 
counterpart, x 2 , is given by 19 

1 d 3 m _ 47rS\ 2 _ 
X2 = 6dlP = ( ^T^ X2 ' 

X2 = ^(^) 3 [<m 4 >-3<m 2 > 2 ]. (2.10) 

In these expressions of \ and X2, we omit the parts which are odd under the time- 
reversal operation, m — >■ — m: In zero external field, such odd parts should vanish 
identically in any finite system in full equilibrium. Note that the above X2, being 
proportional to the minus of the third-harmonic component of the ac susceptibility, 
is sometimes denoted as X3 i n the literature. 



3 Monte Carlo simulation 

We perform Monte Carlo simulations for the Hamiltonian (2.1) based on the stan- 
dard Metropolis method combined with an extended ensemble method. Several 
values of the dimensionless inductance, C, are simulated including £ = 1,3,4,5, 
with the greatest effort in the case of C = 1. Since the bare penetration depth, Ao, 
in those cases is equal to, or less than one lattice spacing, we expect that the effect 
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of screening should manifest itself even for rather small lattices studied here, which 
contain L x L x L sites with L = 3,4, 6, 8, 10. Sample average is taken over 1540 
(L = 3), 1000 (L = 4), 500 (L = 6), 300 (L = 8) and 100 (L = 10) independent 
bond realizations. 

We impose free boundary conditions on all sides of the lattice. If, on the other 
hand, one applies the standard periodic boundary conditions on the link variables 
Ai/s, one has somewhat unphysical results that the magnetization vanishes triv- 
ially even under external fields. In zero field, the periodic boundary conditions 
imposed on the link variables Ai/s also lead to the trivial vanishing of the suscep- 
tibility. In view of such unphysical nature of the periodic boundary conditions, we 
impose free boundary conditions in the following simulations. 

Monte Carlo simulation is performed according to the version of an extended 
ensemble method of Ref.20, where the whole configurations at two neighboring 
temperatures of the same sample are occasionally exchanged. In this method, one 
simulates the sample with a given bond realization at Nt distinct temperatures 
at a time distributed in the range [T min ,T max \. Monte Carlo updating consists of 
the two parts: The first part is the standard local Metropolis updating at each 
temperature, and the second part is an exchange of the whole lattices at two 
neighboring temperatures. 

Since the present spin-glass-like model possesses the link variable, A^, in ad- 
dition to the site variable, Qi, an equilibrium simulation is rather hard even with 
the new efficient algorithm. In the case of L = 8, for example, we prepare 20 
temperature points distributed in the range [0.27J, 0.8J] for a given sample, and 
perform 1.5 x 10 5 exchanges per temperature of the whole lattices combined with 
the same number of standard 'single-spin-flip' Metropolis sweeps. 20 Equilibration 
is checked by monitoring the stability of the results against at least three-times 
longer runs for a subset of samples. 
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As long as one is interested in the gauge-invariant quantities like the chirality 
or the flux, the results would not depend on the particular choice of the gauge. In 
most of our calculation, we choose the gauge where the A^s along the z-direction 
are fixed to be zero. We also use other gauges to take some limited data, including 
the Coulomb gauge and the 'temporal gauge', 18 ' 19 just to make sure that the results 
are really independent of the choice of the gauge. 

We run in parallel two independent replicas with the same bond realization and 
compute an overlap between the chiral variables in the two replicas, 17 

In terms of this chiral overlap, q K , the Binder ratio of the chirality is calculated by 

*» = 5< 3 -R?^>- (3 ' 2) 

Here gcc is normalized so that in the thermodynamic limit it tends to zero above 
the chiral-glass transition temperature, Tec an d tends to unity below Tqg provided 
the ground state is non-degenerate. At the chiral-glass transition point, curves of 
gee against T for different L should intersect asymptotically. 

The chiral-glass susceptibility, which is expected to diverge at the chiral-glass 
transition point, is given by 

Xcg = N p [<q 2 K >}, (3.3) 

The behavior of the chiral short-range order may be seen via the mean magnitude 
of the local chirality, R, defined by 

^ = {^E(<^>]} 1/2 - (3-4) 

Note that for the frustrated noncollinear configurations, R, tends to take a finite 
value whereas for the unfrustrated collinear configurations R tends to vanish. One 
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can also define a reduced chiral-glass susceptibility, xcg, corrected for the short- 
range order effect by dividing xcg by the appropriate power of the magnitude of 
the local chirality, 17 

Xcg = Xcg/^ 4 - (3.5) 

Similarly to the case of the chirality, one can introduce an overlap between the 
flux variables in the two independent replicas, 

?/ = ^E/i 1) /f. (3.6) 

In terms of qf, the Binder ratio of the flux is calculated by 

-A 



2 l [<q}>\- 
The flux-glass susceptibility, xfgj an d its reduced counterpart, xfg; are defined by 



*°=2< 3 -R^# (3.7) 



Xfg = N p [< q 2 f >], (3.8) 

Xfg = Xfg// 4 , (3.9) 

respectively, where the the mean magnitude of the local flux is defined by 

;={^Et</ P 2 >]} 1/2 - (3-io) 



4 Monte Carlo results 

In this section, we present the results of our Monte Carlo simulations. Most ex- 
tensive simulation is made for the inductance C = 1, which corresponds to the 
bare penetration depth, Ao, equal to one lattice spacing. In the subsection [A], we 
present our Monte Carlo results for this inductance, C = 1. The results for other 
inductactances including C = 3,4,5 will be presented later in the subsection [B]. 

12 



[A] I = 1 

In Fig.l, the temperature dependence of the root-mean square of the local- 
chirality amplitude, R, defined by Eq.(3.4), is shown for various lattice sizes. Even 
at lower temperatures R keeps a nonzero value, slightly increasing with decreasing 
temperature, which indicates that the chiral short-range order is developed in the 
temperature range studied here. 

Fig. 2a displays the size and temperature dependence of the Binder ratio of the 
chirality, #cg- The data of gcc for L — 3,4,6,8 all cross at almost the same tem- 
perature T ~ 0.28—0.29, strongly suggesting the occurrence of a finite-temperature 
chiral-glass transition at Tqg — 0.286 ± 0.01 (temperature T is measured in units 
of J). In particular, the data below Tqg show a rather clear fan out. 

The determined value of Tqq is slightly lower than the corresponding chiral- 
glass transition temperature of the pure ±J XY spin glass determined in Ref.17, 
Tqg — 0.32±0.01. Note that the spin-glass model corresponds to the £ — > limit of 
the present model. The observed suppression of Tcg by the screening effect seems 
reasonable, since the screening effect makes the long-ranged interaction between 
vortices (chiralities) short-ranged, which should make the chiral-glass transition 
less favorable. 

The value of the Binder ratio at the transition point, gcQ, * s estimated to be 
g^Q ~ 0.38. The estimated value is considerably smaller than the corresponding 
value of the 3D XY spin glass, 17 g CG ~ 0.72. This large deviation probably 
comes from the difference in the choice of boundary conditions, i.e., free boundary 
conditions in the present simulation and periodic boundary conditions in Ref.17. 
Note that the value of g* is known to depend on the choice of boundary conditions 
even in a given universality class. 

Standard finite-size scaling analysis is made for gcG based on the one-parameter 
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fit of the form, 

gc G = g CG (L l i^ \T-T CG \), (4.1) 

with fixing Tcg = 0.286, where gcc is a scaling function. Then, the chiral 
correlation-length exponent uqg is estimated to be uqg = 1-3 ± 0.2. The cor- 
responding finite-size-scaling plot in given in Fig. 2b. 

The temperature and size dependence of the reduced chiral-glass susceptibility, 
XcG; defined by eq.(3.5), are shown in Fig.3a. Finite-size scaling analysis based on 
the relation, 

Xcg = L 2 '^x CG (L 1/vcG I T - T CG |), (4.2) 

is made with fixing Tcg — 0.286 and vqq = 1.3, yielding the chiral critical-point- 
decay exponent 7/cg — —0.2 ±0.2. The resulting finite-size-scaling plot is displayed 
in Fig.3b. Other exponents can be estimated via the standard scaling relations as 
7 CG ~ 2.9 and (3 CG ~ 0.5. 

The obtained chiral-glass exponents are rather close to the values determined 
previously for the ±J XY spin glass, i.e., the model without screening; vqq = 
1.5 ± 0.3 and rjcc — —0.4 ± 0.2. 17 Therefore, our present result seems consistent 
with the view that the screening effect is actually irrelevant at the 3D chiral-glass 
transition. 

It should be noted here that the determined chiral-glass exponents are also 
not far from the standard spin-glass exponents of the 3D Ising spin glass. 22,23 In 
recent Monte Carlo simulations of the 3D Ising spin glass, however, considerably 
different values were reported for the exponent v, depending on whether v was 
determined from the scaling of the Binder ratio, g, or from the scaling of the 
spin-glass susceptibility, xsg- in Ref.22, for example, the former procedure gave 
an estimate v ~ 2.0 while the latter procedure gave v ~ 1.6, whereas in Ref.23, 
the former gave v ~ 3.0 and the latter gave v ~ 1.5. By contrast, in the present 
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simulation, we did not observe such significant discrepancy between the estimated 
values of vqg'- A common value of vqq ~ 1.3 gave satisfactory fits both for gcc 
and Xcg- At present, we are not sure whether such apparent difference from 
the standard Ising spin glass is simply due to finite-size effects, or is suggesting 
the possibility that the universality class of the chiral-glass transition and that of 
the Ising spin glass are in fact different. Since the chirality can be viewed as a 
"continuous" Ising variable, there exists an obvious similarity between the chiral 
glass and the pure Ising spin glass from a symmetry viewpoint. By contrast, in the 
present model, there exists a local constraint on the possible configurations of the 
chiral variables, 24 which is absent in the Ising spin glass. This may possibly change 
the universality class of the transition. Further studies are required to clarify this 
point. 

Anyway, the occurrence of an equilibrium ordered phase appears to be clear 
from our numerical data, and is in sharp contrast to the vortex-glass (gauge- 
glass) problem where the screening was found to destabilize the equilibrium ordered 
phase. 3 Presumably, such difference comes from the fact that the broken symmetry 
is a discrete Zi symmetry here while it is a continuous U(l) symmetry in Ref.3. 

We also study the behavior of the flux. The temperature dependence of the 
root-mean square of the local-flux amplitude, /, defined by Eq.(3.10), is shown in 
Fig. 4. As can be seen from the figure, the magnitude of the induced local flux is of 
order 0.1 flux quantum per plaquette for this inductance (C = 1). Note that, in the 
small inductance limit C — > 0, / tends to zero, while in the large inductance limit 
C — ► oo, / tends to 1/2 in the ground state of an isolated frustrated plaquette. 

The flux Binder ratio, g-pc, defined by Eq.(3.7), and the reduced flux-glass 
susceptibility, xfg, defined by Eq.(3.9), are shown in Fig.5,and 6, respectively. 
Naively, one expects that the flux should behave in the same way as the chirality, 
since the flux is also a pseudoscalar variable sharing the same symmetry as the 
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chirality. Indeed, as can be seen from Fig. 6, the flux-glass susceptibility, xfgj shows 
a divergent behavior similar to xcg- However, in contrast to the naive expectation, 
clear crossing of the Binder ratio as observed in qqq is not observed in gpc & t least 
in the range of lattice sizes studied here: see Fig. 5. Rather, the ordering tendency 
seems more enhanced in the sense that #fg tends to increase with increasing L 
exhibiting a feature of the ordered phase even above Tqg — 0.286. As the flux 
is an induced quantity generated by the finite inductance effect, we believe this 
behavior to be a finite-size effect. Presumably, for this inductance value, the flux 
hardly reaches its asymptotic critical behavior in rather small lattices studied here. 
Recall that, in the C — ► limit, the flux vanishes trivially while the chiral-glass 
transition itself is most favored. In fact, for larger inductances, we have found that 
the Binder ratios of the flux and of the chirality show more similar behavior as 
expected (see below). 

We also compute the zero-field linear and nonlinear susceptibilities, x an d X2, 
defined by Eqs.(2.9) and (2.10), respectively. As can be seen from Fig. 7, the linear 
susceptibility turns out to be paramagnetic, x > 0, over an entire temperature 
range studied, including in the disordered phase T > Tqg, without a clear anomaly 
at Tqg- in shorter simulations on the same model where the full equilibration is 
not achieved, x tends to get smaller and sometimes becomes negative. 19 These 
results seem consistent with an earlier finding of Dominguez et al who observed a 
paramagnetic x m the field-cooling mode, but a diamagnetic x m the zero-field- 
cooling mode. 18 Meanwhile, the simulation of Ref.18 was performed for a rather 
large inductance, C = 8, where the chiral-glass transition probably did not occur 
in equilibrium (see below). It should be stressed here that the sign of x is i n 
fact a nonuniversal property: Effects not taken into the present model, such as 
the contribution of intragranular supercurrents, could give additional diamagnetic 
contribution in real systems and could easily change the sign of the observed x- 
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By contrast, on general theoretical grounds, the nonlinear susceptibility, X2, 
is expected to show a negative divergence at the transition point where the time- 
reversal symmetry is spontaneously broken in a spatially random manner. 10 Indeed, 
as shown in Fig.8(a), we have observed a behavior fully consistent with this ex- 
pectation. The exponent associated with this negative divergence is estimated via 
a finite-size scaling analysis with assuming Tqg — 0.286 and ucq=1.3, yielding 
72 ~ 4.4 (see Fig. 8b), where \2 ~| T — Tqg |~ 72 - This value of 72 is somewhat 
larger than the chiral-glass susceptibility exponent, 7cg — 2.9. At present, it is 
not entirely clear whether this deviation reflects a true difference in the asymp- 
totic critical behavior. It appears likely that the observed larger value of 72 simply 
comes from a finite-size effect related to the possible nonasymptotic behavior of 
the flux as observed in g F G- 

[B] I = 3, 4 and 5 

So far, the results reported were exclusively for the inductance C = 1. We have 
also made similar, but less extensive simulations for other inductances C = 3,4,5 
in order to study the inductance dependence of the chiral-glass ordering. 

In Figs. 9 and 10, the temperature and inductance dependence of the magnitude 
of the local chirality, R, and that of the local flux, /, are shown for a fixed lattice 
size L = 6. With increasing C, R tends to be suppressed while / tends to be 
enhanced. In Fig.ll(a)-(c), the temperature dependence of the chiral Binder ratio, 
gcG, is shown for C = 3, 4 and 5, respectively. For £ = 3,4, the curves of gcc f° r 
different L appear to cross at a finite temperature. As expected, the chiral-glass 
transition temperature monotonically decreases as C increases. For C = 5, on 
the other hand, no crossing of gcc is observed in the temperature range T > 0.1, 
suggesting that the chiral-glass transition is highly suppressed at this inductance. 
We have also tried similar simulations for even larger inductances, C > 5. However, 
the relaxation becomes extremely slow for these larger inductances, and we can no 
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longer equilibrate the system down to the temperature range of interest within a 
reasonable amount of computation time. 

The tendency that the chiral-glass ordering is suppressed at larger inductances 
can also be seen from other quantities, such as the reduced chiral-glass susceptibil- 
ity, Xcg, the reduced flux-glass susceptibility, xfg, and the nonlinear susceptibility, 
Xi- Thus, we show in Fig.12, 13 and 14 the temperature and inductance depen- 
dence of Xcg, Xfg and \2 for a fixed lattice size L = 6. As can clearly be seen 
in these figures, the chiral-glass ordering is more and more suppressed for larger 
C By contrast, the paramagnetic tendency of the linear susceptibility tends to be 
enhanced for larger C. This is evident from Fig. 15 in which the temperature and 
inductance dependence of the linear susceptibility, x? is shown for a fixed lattice 
size L = 6. 

The obtained phase diagram in the T — C plane is sketched in Fig. 16. There 
appears to be a finite critical value of the inductance, C c , above which there is no 
equilibrium chiral-glass transition. Although it is difficult to precisely locate C c due 
to the extremely slow relaxations we observed at lower temperatures, it appears 
to lie around 5 IS C c IS 7. If this is the case, the value of the inductance chosen by 
Dominguez et al 18 lied in the region of the phase diagram where no equilibrium 
chiral-glass transition took place. Then, a kind of cooperative phenomenon accom- 
panied with a sharp growth of the paramagnetic \i which was observed around 
T ~ 0.4 in Ref.18, might not be related to an equilibrium phase transition, but be 
purely of dynamical origin. This is consistent with our observation in Fig. 15 that 
the paramagnetic tendency is more enhanced for larger C, while the chiral-glass 
ordering itself is suppressed for larger C. 

In Figs. 17(a) and (b), we show the temperature dependence of the flux Binder 
ratio, g-pGi f° r the case of C = 4 and 5, respectively. In the case of C = 4, the curves 
of qfq of L = 4 and L = 6 almost cross, while they do not cross for C = 5. Such 
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behavior is more or less similar to the one observed in the corresponding chiral 
Binder ratio, gcc, for these inductances: see Figs. 11(b) and (c). This observation 
suggests that for these larger inductances the flux behaves in the same way as the 
chirality even in rather small lattices studied here, in contrast to the C = 1 case. 

5 Summary and discussion 

We have shown by extensive Monte Carlo simulations that an equilibrium zero-field 
phase with a spontaneous broken time-reversal symmetry, a chiral-glass phase, is 
possible in certain ceramic superconductors with anisotropic pairing symmetry. 
This phase is truly stable even in the presence of screening. As in spin glasses, the 
nonlinear susceptibility exhibits a negative divergence at the chiral-glass transition 
point. Via a finite-size scaling, static exponents associated with the chiral-glass 
transition are determined. The obtained exponents are rather close to those of 
the 3D Ising spin glass. A rough phase diagram is constructed in the temperature- 
inductance plane. It is found that the chiral-glass transition tends to be suppressed 
for larger inductances, and there appears to be a critical value of the parameter C 
beyond which there is no equilibrium chiral-glass phase. By contrast, the paramag- 
netic tendency of the linear susceptibility (paramagnetic Meissner effect) tends to 
be enhanced for larger inductances. This observation clearly shows that, although 
the paramagnetic Meissner effect could also arise from the frustration effect asso- 
ciated with the anisotropic nature of the superconducting order parameter, it has 
no direct relevance to the chiral-glass phase and the chiral-glass transition. 

Next, on the basis of our findings in the present paper, we wish to discuss 
some requirements for the appropriate ceramic samples where one could expect 
the chiral-glass phase. One important parameter characterizing the sample is the 
dimensionless inductance, C, given by Eq.(2.3). Our present result suggests that an 
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equilibrium chiral-glass state could be realized in the type of samples with smaller 
C, but not for the samples with larger C. If one models a loop as a cylinder of 
radius r and height h, its inductance is given by C = 47i 2 r 2 /h. Putting r ~ l//m, 
h/r ~ 0.01 and J ~ 20K (these values are chosen to mimic the sample used in 
Ref.16), one gets C ~ 10~ 2 . Since this value is considerably smaller than C c , 
an equilibrium chiral-glass phase may well occur in such samples. By contrast, 
if the sample has too large a grain size or too strong Josephson coupling, an 
equilibrium chiral-glass phase will not be realized, or at least largely suppressed. 
Another requirement for the sample is that the grains must be connected via weak 
links into an infinite cluster, not decomposed into finite clusters. Obviously, finite- 
cluster samples cannot exhibit a chiral-glass transition, although the paramagnetic 
Meissner effect is still possible. 6 ' 9 

Once appropriate samples could be prepared, the chiral-glass transition is de- 
tectable in principle via the standard magnetic measurements by looking for a 
negative divergence of \2 or a magnetic ageing phenomenon, as in the case of spin 
glasses. In such magnetic measurements, care has to be taken to keep the intensity 
of applied ac and dc fields weak enough, typically much less than 1G, so that the 
external flux per loop is sufficiently smaller than (p . Recently, a sharp negatively- 
divergent anomaly of \i was reported in a YB2C4O8 ceramic sample by the ac 
method by Matsuura et al, 16 which might be a signal of the chiral-glass transition. 
Ageing was observed in certain ceramic samples, 25 but not in other samples. 9 

As in the case of spin glasses, measurements of dynamic susceptibilities such 
as x'i 00 ) an d x"(w) would also give useful information on the possible chiral-glass 
ordering, particularly when combined with the dynamic scaling analysis. For ex- 
ample, near the chiral-glass transition point, the imaginary part of the linear sus- 
ceptibility, ^'{lo), is expected to satisfy the dynamic scaling relation of the form, 

x"(cv,T,H)^^ cg/zcgucg x"(^^, -Al — ),) (5-1) 
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where t =\ (T — Tqg)/Tcg I an d ^cg is a dynamical chiral- glass exponent. From 
the present calculation, we get the static chiral-glass exponents to be z/qg — 1-3, 
/3 CG ~ 0.5 and j3 CG + 7 CG — 3.4. Although we cannot give a direct numerical 
estimate of the dynamical exponent from the present simulation, one might guess 
that ^cg would take a value around ^cg^cg — 7 — 8 if one would assume the 
analogy between the chiral glass and the Ising spin glass also in the dynamics. 26 

Indeed, a dynamic scaling analysis was made by Leylekian, Ocio and Hammann 
for LSCO ceramic samples. 15 These authors performed both the ac susceptibility 
and the noise measurements, and found an intergranular cooperative transition 
even in zero field at a temperature about 10% below the superconducting transition 
temperature of the grain. Note that the noise measurements enable one to probe 
truly zero-field phenomena where one can be free from the extrinsic pinning effects 
such as the ones envisaged in the so-called critical-state model. 27 It was then found 
that the data of x" satisfied the dynamic scaling of the form (5.1). Here note that 
one is not allowed to invoke the standard vortex-glass scenario to explain such 
intergranular cooperative transition in zero field, since in the standard vortex- 
glass picture frustration is possible only under finite external fields. By contrast, 
the experiment seems consistent with the chiral-glass picture. 

Meanwhile, when the intragranular superconducting transition and the inter- 
granular transition take place at mutually close temperatures as in Ref.15, the 
Josephson coupling, J, which has been assumed to be temperature independent 
in the present model, is actually strongly temperature dependent in the transition 
region. In such a case, care has to be taken in analyzing the experimental data, 
since such temperature dependence of J might modify the apparent exponent value 
from the true asymptotic value to some effective value. In fact, the dynamical ex- 
ponent zv ~ 30 determined by Leylekian et al were different from the standard 
spin-glass value, which might be due to the proximity effect of the intragranular 
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superconducting transition. 15 If, on the other hand, the temperature dependence 
of J was taken into account in the fit, a more realistic value zv ~ 10 — 15 was ob- 
tained in Ref.15. By contrast, when the intergranular chiral-glass transition takes 
place at a temperature much lower than the intragranular superconducting tran- 
sition temperature as in Ref.16, the Josephson coupling can safely be regarded as 
temperature independent as assumed in the present model. Anyway, it is desirable 
to get a direct numerical estimate of the dynamical chiral-glass exponent, zcg 5 to 
be compared with experiments. We are now planning to perform a simulation to 
get an independent numerical estimate the dynamical exponent. 

It may also be possible to detect a spontaneously induced flux in the chiral- 
glass state by muon spin relaxation or electron holography in zero external field. 
As in the noise measurements, these measurements can be made in zero external 
field, and has an advantage of being free from the pinning effects of extrinsic 
origin. Here it is essential to make such measurements for ceramic samples with 
sufficiently many weak links, not for single crystals, simply because the chiral-glass 
phase is expected only in the former type of samples. By contrast, the kind of 
time-reversal-symmetry breaking proposed by Wen, Wilczek and Zee is associated 
with the time-reversal-symmetry of the bulk superconducting order parameter and 
should occur even in single crystals. 28 

We could estimate an order of the induced flux below Tqg from the results of 
our present simulation. In the case of C = 1, for example, the flux intensity can 
be estimated from the calculated [< qj >] to be about O.O20 O at 20% below Tqg- 
For a sample with a typical grain diameter about 1/im, this corresponds to the 
flux intensity equal to 0.4G, which seems well within the sensitivity of the yUsR 
measurements. For a sample with a grain diameter about 10//m, the flux intensity 
would be reduced to 4mG. If the dimensionless inductance, C, is significantly 
smaller than unity, the flux intensity would become much smaller, and eventually 
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vanishes in the C — > limit. 

In the chiral-glass state, the U(l) gauge symmetry will not be broken, even 
randomly, in the strict sense. This means that the phase of the condensate, 9, 
remains disordered in the chiral-glass state at least on sufficiently long length and 
time scales. Thus, the chiral-glass state should not be a true superconductor, with a 
small but nonvanishing linear resistance even below Tqg- This property has been 
established in the C — » limit where the screening effect can be neglected. 10,17 
Although we have not measured in the present simulation the quantity directly 
relevant to the U(l) gauge-symmetry breaking, the screening effect makes the 
interaction between vortices short-ranged and makes the the U(l) gauge-symmetry 
breaking transition even more unlikely. 

Small but finite linear resistivity, p^, in the chiral-glass state arises due to the 
slow motion of thermally-activated integer vortex lines (unbound vortex loops). 
Free motion of integer vortex lines is still possible in the chiral-glass state where 
chiralities (half-vortices) sitting at frustrated plaquettes are frozen. A schematic 
picture showing such free motion of integer vortex-line excitations in the back- 
ground of a frozen pattern of chiralities is given in Fig. 18. One can see that free 
motion of integer-vortex lines of either sign is possible without seriously destroy- 
ing the freezing pattern of chiralities in the background. In order to destroy the 
chiral-glass ordering in the background, a chiral domain-wall-type excitation is 
necessary, which would be responsible for the chiral-glass transition at T = Tcg- 
On decreasing the temperature across Tqg, a sharp drop of the resistivity will be 
caused by such chiral domain-wall excitations, but the resistivity will stay finite 
even below Tcg due to the wandering vortex-line excitations. 

We try to get a very rough order estimate of pl at the chiral-glass transition 
point based on a flux-creep model. 29 Within this model, the linear I — V relation 
with finite pi is expected below a characteristic current-density scale, j c , given by 
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j c ~ A;bT/(0o^ 2| C 2 )) where d is a typical grain size and £ is the phase correlation 
length (or the 'spin' correlation length) in units of d. We estimate £ at T = Tqq 
from the Monte Carlo data of the 3D XY spin glass 30 ' 17 to be £ ~ 10 lattice 
spacings. Here note that the spin-glass correlation length does not diverge at 
the chiral-glass transition point. 17 Then, for a typical ceramic sample, we put 
d ~ 1/im, T = T CG ~ 30K, to get j c ~ 2 x 10 3 A/m 2 . The flux-creep model 
also yields pl ~ ^^/{ksTTof), where To is an inverse 'attempt frequency' of 
the intergranular vortex motion, and f is the phase or 'spin' correlation time in 
units of To. Again, from the Monte Carlo data, 30 ' 17 we estimate f at T = Tqg 
to be f ~ 5 x 10 4 Monte Carlo time steps. Precise value of our time unit, r , is 
largely unknown, but it should be much longer than the atomic time scale since 
the vortex motion of interest here is the one over grains. If we put r ~ 10~ 9 sec, 
for example, we have p^ ~ 0.2/ifi-cm at T = Tqg, while for tq ~ 10~ 5 sec, we have 
p L ~ 0.2 x 10~ 4 /if2-cm. These values, though small, may be within the reach of 
careful experimental measurements. 

All simulations presented in this paper were done in zero external field. A 
chiral-glass phase and a chiral-glass transition are associated with a spontaneously 
breaking of time-reversal symmetry, and in that sense, can be regarded as a zero- 
field phenomenon. Still, it should be emphasized here that the fate of the chiral- 
glass phase and the chiral-glass transition in an external field is not necessarily 
trivial and is of great interest. Clearly, under external magnetic fields, the system 
no longer possesses a global time-reversal symmetry. Therefore, there cannot be 
a chiral-glass transition associated with a spontaneous breaking of a global time- 
reversal symmetry. Nevertheless, an interesting possibility emerges if the chiral- 
glass transition in zero field accompanies the replica symmetry breaking 31 of the 
chirality. In such a case, an equilibrium phase and the associated thermodynamic 
transition should persist even under finite fields and are characterized by the chiral 
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replica symmetry breaking. Then, the transition line in the H — T plane might look 
like the so-called AT-line familiar in spin glasses. 32 It is very interesting to relate 
such chiral-AT line to the AT-like line often observed experimentally in ceramic 
high-T c superconductors. 12 

We finally note that the chiral-glass state can be realized not only in high-T c 
superconductors, but also in other superconductors with nontrivial pairing sym- 
metry, such as in heavy fermion superconductors, or possibly, in some organic 
superconductors. It would be interesting to experimentally search for this novel 
phase in these materials, since it is a new state of matter realized only in anisotropic 
superconductors with unconventional pairing symmetry. 

The numerical calculation was performed on the FACOM VPP500 at the su- 
percomputer center, Institute of Solid State Physics, University of Tokyo. A part 
of this work was made when one of the authors (M.S. Li) was in Kyoto Institute 
of Technology. M.S. Li thanks the Japan Society for Promotion of Science for the 
award of a fellowship. He was also supported in part by the Polish KBN grant. 



25 



References 

1 G. Blatter, M.V. Feigel'man, V.B. Geshkenbein, A.I. Larkin and V.M. Vinokur, 

Rev. Mod. Phys. 66 (1994) 1125. 

2 M.P.A. Fisher, Phys. Rev. Lett. 62 (1989) 1415. 

3 H.S. Bokil and A.P Young, Phys. Rev. Lett. 74, (1995) 3021; C. Wengel and 

A.P Young, Phys. Rev. B54 (1996) R6869. 

4 See, for example, A. Mathai, Y. Gim, R.C. Black, A. Amar and F.C. Wellstood, 

Phys. Rev. Lett. 74 (1995) 4523; D.J. van Harlingen, Rev. Mod. Phys. 67, 
515 (1995). 

5 F.V. Kusmartsev, Phys. Rev. Lett. 69 (1992) 2268; J. of Superconductivity 5 

(1992) 463. 

6 M. Sigrist and T.M. Rice, J. Phys. Soc. Jpn. 61 (1992) 4283; Rev. Mod. Phys. 

67 (1995) 503. 

7 P. Svedlindh, K. Niskanen, P. Nordblad, L. Lundgren, B. Lonnberg and T. 

Lundstrdm, Physica C162-164 (1989) 1365. 

8 W. Braunisch, N. Knauf, G. Bauer, A. Kock, A. Becker, B. Freitag, A. Griitz, 

V. Kataev, S. Neuhausen, B. Roden, D. Khomskii, D. Wohlleben, J. Bock 
and E. Preisler, Phys. Rev. B48 (1993) 4030. 

9 J. Magnusson, J.-O. Andersson, M. Bjornander, P. Nordblad and P. Svedlindh, 

Phys. Rev. B51 (1995) 12776; J. Magnusson, M. Bjornander, L. Pust, P. 
Svedlindh, P. Nordblad and T. Lundstrom, Phys. Rev. B52 (1995) 7675. 

10 H. Kawamura, J. Phys. Soc. Jpn. 64 (1995) 711. 



26 



11 Note that the same state was also called 'orbital-glass phase' in Ref.10. Since 

the term 'orbital-glass' has often been used in the literature simply to de- 
scribe the state exhibiting the paramagnetic Meissner effect, however, we 
donot use this terminology here to represent the thermodynamic phase of 
our interest, but rather, use the term 'chiral-glass', following the XY spin- 
glass terminology. In the chiral-glass phase, the susceptibility could be either 
paramagnetic or diamagnetic. 

12 K.A. Muller, M. Takashige and J.G. Bednorz, Phys. Rev. Lett. 58 (1987) 

1143. 

13 Z. Koziol, Physica C159 (1989) 281. 

14 K. Park, J.J. Kim and J.C. Park, Solid State Comm. 71 (1989) 743. 

15 L. Leylekian, M. Ocio and J. Hammann, Physica C185-189 (1991) 2243; 

Physica B194-196 (1994) 1865. 

16 M. Matsuura, M. Kawachi, K. Miyoshi, M. Hagiwara and K. Koyama, J. Phys. 

Soc. Jpn. 64 (1995) 4540. 

17 H. Kawamura, Phys. Rev. B51 (1995) 12398;. J. Phys. Soc. Jpn. 61 (1992) 

3062; H. Kawamura and M. Tanemura, J. Phys. Soc. Jpn. 60 (1991) 608. 

18 D. Dommguez, E.A. Jagla and C.A. Balseiro, Phys. Rev. Lett. 72 (1994) 

2773. 

19 H. Kawamura and M.S. Li, Phys. Rev. B53 (1996) 619. 

20 K. Hukushima and K. Nemoto, J. Phys. Soc. Jpn. 65 (1996) 1604. 

21 H. Kawamura and M.S. Li, to appear in Phys. Rev. Lett. (1997). 



27 



22 N. Kawashima and A.P. Young, Phys. Rev. B53 (1996) 484; A.P. Young and 

N. Kawashima, Int. J. Mod. Phys. C7 (1996) 327. 

23 K. Hukushima, H. Takayama and K. Nemoto, Int. J. Mod. Phys. C7 (1996) 

337. 

24 J. Villain, J. Phys. Cll (1978) 745. 

25 C. Rossel, Y. Maeno and I. Morgenstern, Phys. Rev. Lett. 62 (1989) 681. 

26 A.T. Ogielsky, Phys. Rev. B32 (1985) 7384. 

27 C.P Bean, Rev. Mod. Phys. 36 (1964) 31. 

28 X.G. Wen, F. Wilczek and A. Zee, Phys. Rev. B39 (1989) 11413. 

29 PW. Anderson and Y.B. Kim, Rev. Mod. Phys. 36 (1964) 39. 

30 S. Jain and A.P. Young: J. Phys. C19 (1986) 3913. 

31 G. Parisi, J. Phys. A13 (1980) 1101. 

32 J.R.L. de Almeida and D.J. Thouless, J. Phys. All (1978) 983. 



2S 



Figure captions 

Fig. 1 The temperature and size dependence of the root-mean square of the local- 
chirality amplitude, R, defined by Eq.(3.4), for £ — 1. 

Fig. 2 (a) The temperature and size dependence of the Binder ratio of the chirality, 
gcc f° r C — 1. Inset is a magnified view around the transition temperature 
Tqg — 0.286. (b) Finite-size scaling plot of gcc with Tqg = 0.286 and 
z/cg = 1-3. 

Fig. 3 (a) The temperature and size dependence of the reduced chiral-glass suscep- 
tibility, XcG> f° r £ — 1- (b) Finite-size scaling plot of x CG with T CG = 0.286, 
z/cg = 1-3 and r] CG = -0.2. 

Fig. 4 The temperature and size dependence of the root-mean square of the local- 
flux amplitude, /, defined by Eq.(3.10), for C = 1. 

Fig. 5 The temperature and size dependence of the Binder ratio of the flux, gpQ, 
for £ = 1. 

Fig. 6 The temperature and size dependence of the reduced flux-glass susceptibility, 
Xfg, for I = 1. 

Fig. 7 The temperature and size dependence of the zero-field linear susceptibility, 
X, for C — 1. An arrow in the figure represents the location of the chiral-glass 
transition point. 

Fig. 8 (a) The temperature and size dependence of the zero-field nonlinear suscep- 
tibility, X2, for £ — 1. (b) Finite-size scaling plot of X2 with Tqg = 0.286, 
z/cg — 1-3 and 72 = 4.4, where X2 ~ (T — Tcg) 72 - An arrow in the figure 
represents the location of the chiral-glass transition point. 
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Fig. 9 The temperature and inductance dependence of the root-mean square of the 
local-chirality amplitude, k, defined by Eq.(3.4), for a fixed lattice size L = 6. 

Fig. 10 The temperature and inductance dependence of the root-mean square of the 
local-flux amplitude, /, defined by Eq.(3.10), for a fixed lattice size L = 6. 

Fig. 11 The temperature and size dependence of the Binder ratio of the chirality, 
gee, f° r ( a ) £ — 3, (b) C = 4 and (c) C = 5. 

Fig. 12 The temperature and inductance dependence of the reduced chiral-glass 
susceptibility, Xcgj f° r a fixed lattice size L = 6. 

Fig. 13 The temperature and inductance dependence of the reduced flux-glass sus- 
ceptibility, xfg; f° r a fixed lattice size L = 6. 

Fig. 14 The temperature and inductance dependence of the zero-field nonlinear sus- 
ceptibility, X2, for a fixed lattice size L = 6. 

Fig. 15 The temperature and inductance dependence of the zero-field linear suscep- 
tibility, X) f° r a fixed lattice size L = 6. 

Fig. 16 A phase diagram in the T-C plane. Renormalized inductance C is defined 
by Eq.(2.3). The data point at C = is taken from Ref.17. 

Fig. 17 The temperature and size dependence of the Binder ratio of the flux, gpQ, 
for (a) C = 4 and (b) C = 5. 

Fig. 18 Two-dimensional segment of the lattice showing thermally-activated integer 
vortex lines with vorticity ±1, wandering in the background of a frozen pat- 
tern of chiralities in the chiral-glass state. Plus (+) and minus (— ) chirality 
can be viewed as half-vortices with vorticity ±1/2 sitting at frustrated pla- 
quettes, while unfrustrated plaquettes are frozen into the zero-chirality (0) 
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state. If one looks at a given frustrated plaquette frozen into the + chiral- 
ity (or vorticity +1/2) state, for example, its vorticity occasionally becomes 
+3/2 or —1/2 when the thermally-activated integer vortex line of either sign, 
+ 1 or -1, passes this plaquette. Still, the long-time average of the vorticity 
at this plaquette is equal to +1/2, showing that the free motion of integer 
vortex lines is compatible with the long-range chiral-glass order. 
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